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ABSTRACT 



NASA’s Plasma Motor Generator (PMG) tethered satellite mission was launched 
in June 1993 to verify the ability of hollow cathode plasma sources to couple electric 
currents from an electrodynamic tether into the ambient ionospheric plasma. This large- 
scale coupling process resulted in turbulent plasma signatures associated with the orbiting 
plasma generator, which propagated over great distances along the earth’s geomagnetic 
field lines. VHF radars in Hilo, Hawaii and Jicamarca, Peru recorded observations of 
these field-aligned disturbances as part of the experiment. Based on analysis of these 
radar observations and tracking data of PMG’s orbit, the effective propagation velocity of 
these traveling plasma waveforms was calculated to be of the order of 1000 meters per 
second. Detection of these disturbances, associated with PMG’s passage overhead, 
supports the existence of a phantom current loop allowing current flow along the 
magnetic field lines of the earth and into the lower ionosphere from either end of an 
electrodynamic tether. 
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I. INTRODUCTION 



At approximately 13:27 GMT on June 26, 1993, an Air Force Delta-II rocket 
roared off the launch pad at Cape Canaveral, its third stage carrying a Global Positioning 
Satellite (GPS) as the primary payload. Attached to the second stage of the Delta-II was 
NASA's Plasma Motor Generator (PMG) tethered satellite mission. Following third stage 
separation and second stage fuel depletion bum, the PMG system was left in an elliptical 
orbit, approximately 207 x 922 kilometers, at a 27.5 degree inclination. At approximately 
2000 seconds mission elapsed time, a Far End Package (FEP) was spring-ejected from the 
Delta-II second stage, at a rate of two to three meters per second, trailing a 500-meter 
conducting wire from the PMG tether deployer assembly. Within minutes, the two-body 
tethered satellite system stabilized in a gravity-gradient configuration with the FEP above 
the Delta-II second stage. 

The primary objective of the PMG experiment was to verify the ability of hollow 
cathode plasma sources to couple electric currents from either end of a long wire moving 
through the ambient Low Earth Orbit (LEO) ionospheric plasma. By providing a 
low-impedance electrical connection with the ambient environment, while moving at 
orbital velocity through the earth’s magnetic field, significant currents can be induced in a 
single-conductor tether wire, providing the capability for direct conversion of orbital 
energy to extractable electrical current (generator principle). Alternatively, orbital 
"steering" forces can be induced through the tether by driving current in the opposite 
direction with applied voltage (motor principle). The PMG experiment utilized these 
physical principles to investigate the overall practicality of reciprocal conversion of orbital 
energy and electricity, and also, to evaluate current technology for transferring charge 
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charge between a spacecraft and the LEO environment, utilizing hollow cathode plasma 
contactors. 

The PMG electrodynamic tether orbited within the tenuous, electrically charged 
layer of the atmosphere known as the ionosphere, providing the opportunity to increase 
understanding of physical processes in the near-Earth space environment, while at the 
same time exploring the mechanical dynamics of tethered systems. A key mission 
objective was to investigate large-scale coupling processes and turbulent signatures 
associated with an orbiting artificial plasma generator and the ambient ionospheric 
plasma. The approach to accomplish this study utilized sensitive VHF meter-wave 
ionospheric radars to sense field-linked plasma perturbations at large distances from the 
orbiting plasma source. Near Lima, Peru, the powerful fifty megahertz radar at Jicamarca 
Radio Observatory was aligned to observe field-aligned backscatter from the near 
horizontal magnetic field lines. Two less sensitive fifty megahertz wind-profilers, located 
at Piura, Peru and Pohnpei, Micronesia, were configured for ionospheric measurements 
along the geomagnetic field lines. Additionally, NASA/JSC deployed a transportable 
radar system to Hilo, Hawaii, developed for remote field study of ionospheric 
modification effects, spacecraft interactions with the lower ionospheric environment, and 
reentering orbital debris. 

Electromagnetic scattering from strong ionospheric turbulence generated during 
operation of the PMG experiment was observed from several of the VHF radar sites by 
radar probing of the geomagnetic flux tubes that mapped from the orbiter through the 
radar's fixed field-of-view. On several passes, field-aligned backscatter was observed 
from plasma structure propagating along the lines at distances up to several hundred 
kilometers from the PMG plasma source. The purpose of this thesis is to present these 
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observations, and develop the implications for the generation of field-aligned 
disturbances by the tether. 
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II. BACKGROUND ON TETHERS 



"Come, let us build ourselves a city, with a tower that reaches to the heavens..." 

— The Tower of Babel, from the Book of Genesis, Chapter 1 1 



A. HISTORY 

The earliest concepts for the application of tethers in space can be traced back to 
the last century. In 1895 Kostantin Tsiolkovsky, the Russian schoolteacher and pioneer of 
astronautics, reflected on a way to create an environment devoid of terrestrial gravity 
through the construction of an equatorial tower pointing into space, reaching beyond 
geostationary altitude. He postulated that upon ascending such a "Heavenly Tower," the 
effects of gravity would gradually diminish until they disappeared completely. Although 
he did not realize it, Tsiolkovsky had conceived the first geostationary 
satellite.(Tsiolkovsky, 1895) 

Sixty-five years later, at the dawn of the space age, the Leningrad engineer Yu 
Artsutanov conceived an alternate concept, consisting of a platform anchored in space, 
pointing towards the earth, with a cable deployed down from the satellite to connect it to 
the earth's surface at the equator. A second cable deployed upward beyond the satellite 
would carry a ballast to maintain the system's center of gravity in the geostationary orbit. 
Artsutanov visualized his "Heavenly Funicular" as a means of supporting transportation 
into space, and thus was bom the idea of a space elevator.( Artsutanov, 1 960) 

The early 1960's saw a flurry of tether-related proposals, including concepts such 
as rolling satellites, lunar and asteroid tethers, space elevators, a space "necklace" around 
the earth, and tether rescues of stranded astronauts. Several early concepts involving 
gravity-gradient stabilized tethers are shown in Figure 1. Initial application experiments 



5 



for tethered vehicle operations in space began in late 1966 as part of the Gemini program. 
The effects of firing translational thrusters and the effects of gravity gradient stabilization 
on the motion of the tethered system were examined using the Gemini XI and XII 
spacecraft with the spent stage of their Atlas-Agena D launch vehicle.(Lang and Nolting, 
1967) Further development of the gravity-gradient concept came in 1969 when A.R. 
Collar and J.W. Flower proposed a very long tether connecting two satellites - one 
beyond geostationary altitude and one in low earth orbit, such that the center of gravity 
was located at geostationary distance. (Collar and Flower, 1969) Additionally, the idea for 
a space propulsion system was introduced by R.D. Moore, consisting of a conducting 
wire with plasma contactors at either end. His "geomagnetic thruster" was a true 
forerunner of present day electrodynamic tethers. (Moore, 1966) 

It was not until the early 1970's, however, that coordinated efforts were made to 
investigate the complex dynamics of long tethers as an innovative approach to space 
transportation. Mario Grossi's proposal of a shuttle-borne tether for use as an ULF/ELF 
antenna and radio physics facility was one of the first practical applications study for the 
use of tethers in space. (Grossi, 1973) Through his involvement in this program, 
beginning in 1974, Giuseppe Colombo quickly recognized the far-reaching potential of 
the tether in space transportation. (Colombo et al., 1974) Professor Colombo, who is 
considered to be the modern pioneer of tethers applications, conceived, introduced or 
analyzed virtually all of the first space propulsion schemes involving tethers. It was these 
new methods of generating propulsion through momentum exchange between spacecraft, 
and in generating electric power and thrust through electrodynamic interactions with the 
space environment, that established tethers in space as a revolutionary operational 
capability today. Tether applications in space exist in a wide range of fields, including 
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aerodynamics, gravity control, electrodynamics, transportation, science studies, and 
planetary exploration. 

B. GRAVITY-GRADIENT STABILIZATION 

There are two basic configurations that provide controlled acceleration fields for 

an orbiting tether system. A gravity-gradient stabilized tether rotates once per orbit in an 

inertial coordinate frame, while a rotating tether configuration consists of multiple 

rotations per orbit. An understanding of the basic physical principles behind gravity- 

gradient forces is fundamental when studying general tether applications. Consider a 

single satellite in a circular orbit about the earth. It is subject to a gravitational force, 

_ - GMm 

F ’=— 

where G is the gravitational constant, M is the mass of the earth, m is the mass of the 
satellite, and r is the radial distance from the center of the earth to the satellite. This force 
must equal the centripetal force necessary to cause circular motion, 

F c = mrcj 2 

where to is the satellite's orbital angular velocity. The balance between the force of 
gravity, which pulls the satellite towards earth, and the outward centrifugal force, keeps 
the satellite in orbit. To remain in a circular orbit, the net vertical force acting on the 
satellite must be zero. Thus, by equating the two forces, the constant angular velocity of 

the satellite as it orbits the earth can be solved for, 

2 GM 
or = — 

r 

Now consider two masses connected by a tether in a circular orbit such that the 
system is oriented vertically. Gravity, centrifugal force, and atmospheric drag all vary 
with altitude, thus each body in a tethered system is subject to separate influences. A 
simple "dumbbell'' configuration, shown in Figure 2, will be in an equilibrium state 
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orbiting the earth as long as the tether remains aligned in the vertical and all parts of the 
system have the same orbital angular velocity. For this orientation, the sum of the 
gravitational forces acting over the system must be equal and opposite to the sum of the 
centrifugal forces acting on the system. Neglecting the tether’s weight, the summation of 
force becomes, 

GMm. GMm , , , 

2 I j — = m \ r t 60 +m 2 r 2<° 

r \ h 

From this equation the tethered system's orbital angular velocity can be computed. The 
gravitational and centrifugal forces acting on the tethered system are equal and balanced 
at only one place - the system's center of gravity. This equilibrium point is in a zero-g 
condition and experiences no net force in the radial direction. It is in free fall as it orbits 
the earth, but the two end masses are not. They are constrained by the tether to orbit with 
the same angular velocity as the center of gravity. The value of the orbital radius for the 
center of gravity is given by, 

, GM 

r » = v 

For short tethers, this value will be nearly equal to the position of the center of mass of 
the system. However, since gravitational acceleration changes non-linearly with distance 
from the center of the earth, the center of gravity of the tethered system does not 
necessarily coincide with the center of mass. The separation is not dramatic for systems 
using short tether lengths, but can be significant for long tethered systems.(Amold, 1987) 
Without the connecting tether, the higher mass would tend to move at a slower 
orbital velocity, while the lower mass would tend to move at a greater orbital velocity. 
Thus, the presence of the tether speeds up the higher mass and slows down the lower 
mass. This causes the higher mass to be subject to a larger centrifugal than gravitational 
force, and the lower mass to be subject to a larger gravitational than centrifugal force. 
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With a tether in place, the net effect of these unbalanced forces is to create tension in the 
tether. The resulting upward acceleration of the higher mass and downward acceleration 
of the lower mass gives rise to a force couple applied to the system, which forces it into a 
vertical orientation. Should a disturbing torque act on the tethered system and displace it 
from the vertical position, this force couple produces restoring forces at each mass, 
illustrated in Figure 3, which act to return the system to a vertical orientation. These 
restoring torques act in both in-plane and out-of-plane directions to tend to keep the 
system in its equilibrium state. 

Although this vertical orientation is a stable one, there are weak but persistent 
forces which cause the tethered system to oscillate about the vertical. These forces 
include atmospheric drag, solar heating effects, electrodynamic forces (for conducting 
tethers), solar pressure torque, and geomagnetic torques. For short tether lengths, where 
the gravity-gradient is relatively constant, the libration frequencies are independent of 
tether length due to the fact that both the displacement and the restoring forces on the 
system vary linearly with tether length. This is not the case for very long tethers, 
however, since the gravity-gradient can vary significantly over the length of the tether. 
The natural restoring torque for out-of-plane disturbances tends to be stronger than that 
for in-plane disturbances, thus out-of-plane librations have a higher frequency than in- 
plane librations. The natural frequency for in-plane oscillations is given by, 
f =2(0 

For out-of-plane oscillations, 

/ = V3 co 

where co is the orbital angular velocity of the center of mass. 

Optimum design considerations for tethers used in space call for minimizing the 
tether's mass while maintaining its required strength. This is accomplished by using 
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materials such as aluminum wire and Teflon insulation which have high strength to 
weight ratios. Additionally, aluminum’s conductivity per mass is about twice that of 
copper, and Teflon provides good resistance to atomic oxygen erosion. The tether cannot 
be too thin because of the significant risk of being cut by a micrometeorite impact. 
Collisions with orbital debris pose a threat to any deployed tether, since any such 
encounter with particles or objects of more than a few grams will almost certainly sever 
it. The tether must be constructed to survive most impacts so as to have a reasonable 
expected lifetime. Tethers designed with a constant cross-sectional area have a limited 
length, however, tethers tapered to maintain a constant stress per unit cross sectional area 
can theoretically have unlimited length. Thus the optimum design for high tether tensions 
would be an exponentially tapered tether with a maximum area at the center of gravity 
and minima at the end masses. Although there is no limit to the length that can be 
achieved, the payload that can be supported at the ends keeps getting smaller.(NASA 
Report, 1989) 

C. ELECTRODYNAMIC TETHERS 

One of the most rewarding applications of tethers to enhance man's use of space 
may be the electrodynamic tether, a satellite system in which two separate payloads are 
connected with an insulated conductive wire. Electrodynamic tethers make use of the 
interactions between a moving conductor, the earth's magnetic field, and the ambient 
plasma in the ionosphere to allow propulsion and power generation. Such a system could 
be used with solar arrays to offset drag in LEO, replace batteries for power storage, or 
provide propulsion for orbital maneuvers. 

Consider the ionospheric plasma, co-rotating with the Earth, relative to which a 
tether is moving with orbital velocity v. For most tether applications v is Eastward, while 
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the earth's magnetic field B, acts Northward. This produces an induced electric field 
acting upward, causing a sizable potential difference across the conducting tether due to 
its motion through the geomagnetic field. By providing contactors (devices capable of 
passing sufficient current to and from the plasma with little voltage drop) at both ends of 
the tether, electrical contact with the Earth's LEO plasma environment is made, 
establishing a current loop through the tether, external plasma, and lower ionosphere. In 
the undriven state, electrons from the ambient plasma are collected by the upper end of 
the tether and emitted from the lower end, creating a net positive charge cloud at the top 
end and a net negative charge cloud at the bottom end. In a quasi-static model, excess 
free charges are constrained to move along the geomagnetic field lines intercepted by the 
ends of the tether until they reach the E region of the lower atmosphere where their are 
sufficient collisions with neutral particles to allow the electrons to migrate across the field 
lines and complete the circuit. Taking satellite motion into account, however, one finds 
that by the time the circuit "closes", the satellite has traveled many kilometers. Thus the 
true dynamic problem appears to be a complex one. Figures 4 and 5 illustrate how 
current flowing along the wire due to the potential difference is coupled into and through 
the ambient ionospheric plasma.(Bainum et al., 1986) 

There are three basic plasma contactor configurations to provide electrical closure 
to the space environment. The first consists of passive large area conductors at both ends 
of the tether. The upper conductor collects electrons, while the lower plasma contactor 
utilizes its large surface area to collect ions. A large conducting balloon could be used at 
the upper end in this configuration, with the conductive surface of the spacecraft serving 
at the lower end. 

An alternate method of current coupling consists of a passive large area conductor 
at the positive end of the tether, and an electron gun at the negative end. In this case, the 
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ejection of a negative electron at a high energy is the equivalent to collecting a positive 
ion current. While this configuration allows higher currents to be achieved, it also 
requires an on-board electrical power supply to drive the electron gun. One general 
problem with electron guns, is that they tend to become current limited, typically by space 
charge effects. 

The final configuration consists of a plasma-generating hollow cathode assembly 
at each tether end. Instead of relying on a passive and physically large conducting surface 
to collect currents, hollow cathodes generate an expanding cloud of highly conductive 
plasma. This plasma provides the necessary thermal electron density to carry the full 
tether current in either direction from a tether end, until it is merged into the ambient 
ionospheric plasma currents. Although this configuration also requires an on-board 
power supply, the bipolar nature of the emitted plasma offers greater flexibility and has 
been demonstrated to work. (Olsen, 1981) 

Hollow cathode systems are considered to be the most desirable contactor 
configuration. Their current reversibility allows the tether to function alternately as either 
a generator or a motor. Additionally, they operate at lower voltages while still producing 
the desired electrical power levels. They are safer for the spacecraft system because they 
establish a vehicle ground reference potential with respect to the local plasma. Finally, 
with hollow cathode contactors, shorter tethers may be used for required current levels, 
greatly reducing the size requirements for stabilizing end masses, and simplifying tether 
deployment and dynamics. (NASA Report, 1989) 

Tethered systems in LEO generate electrical power at a rate equal to the loss in 
spacecraft orbital energy due to drag forces. Consider a vertical, gravity-gradient 
stabilized, insulated, conducting tether, terminated at both ends by hollow cathode plasma 
contactors. As this system orbits the earth, it cuts across the geomagnetic field lines at 
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orbital velocities of about eight kilometers per second. A Faraday electromotive potential 
(emf) is induced in an element of the length of the tether, given by the equation: 

d<t>= (vxB)»Il 

where dO is the induced emf across the tether element length, v is the tether velocity 
relative to the geomagnetic field, B is the magnetic field strength, and dl is a differential 
element of the tether length pointing in the direction of current flow. This emf creates a 
potential difference across the tether by making the upper end of the tether positive with 
respect to the lower end, allowing current to flow through the wire. By using variable 
resistance or impedance matching techniques, the current passing through the tether can 
be controlled, providing power for on-board systems, as shown in Figure 6. 

When the current generated by the induced emf is allowed to flow, a force is 
exerted on the current, and thus on the tether, by the geomagnetic field given by: 

dF= iJlxB 

where dF is the force exerted on an element of the tether by the magnetic field, I is the 
tether current, dl is a differential element of tether length pointing in the direction of 
positive current flow, and B is the magnetic field strength. In LEO, where the relative 
velocity between the orbiting tether and the rotating geomagnetic field is large, this force 
is a drag on the tether. A direct result is that when electrical power is generated by the 
system for on-board use, it is generated at the expense of orbital energy. If the tether is to 
maintain its altitude, this electrodynamic braking must be compensated for by rockets or 
other propulsion means. 

Alternately, when a current from an on-board power supply is fed into the tether 
against the emf induced by the geomagnetic field, the direction of the force is reversed 
and it becomes a propulsive force acting on the tethered system. Figure 7 illustrates this 
principle, which allows the tether can be used as an electrodynamic thruster, with 
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applications of orbital maneuvering and drag compensation. The price to be paid in this 
instance is that the propulsive force is generated at the expense of on-board electrical 
power. 

Electrodynamic tethers also have applications in power storage. Using solar 
arrays, current could be fed into a tether during periods of array illumination, providing a 
propulsive force to boost orbital altitude. Then, during periods of darkness, orbital 
energy could be traded for useful DC electrical power as the geomagnetic field induced 
voltages across the tether, inducing current flow in the opposite direction. This reversible 
energy storage system would complement systems employing the charging and 
discharging of batteries. By itself, however, it is not mass-competitive with conventional 
battery systems. (NASA Report, 1989) 

Thus the electrodynamic tether system presents a new form of continuous power 
or thrust generation in space, which can in theory produce either an efficient power 
supply of a few to hundreds of kilowatts of usable electricity, or provide orbital 
maneuvering forces for the satellite. There are effects due to orbit altitude and inclination 
which act to limit the effective operating range of the electrodynamic tether system to 
orbit altitudes of less than 2000 kilometers and orbit inclinations below sixty degrees. 
The major cause for the altitude effect is the decrease in geomagnetic field strength which 
goes as the inverse of the third power of the orbit radius. A smaller cause is the orbit 
velocity of the system which decreases as the inverse square root of the orbit radius. As 
the orbit inclination is increased, the angle between the orbit velocity vector and the 
magnetic field vector is decreased until at a polar orbit velocity is parallel to the field. At 
this point the vector cross product goes to zero and no voltage is induced in the tether. 
Thus equatorial and low inclination orbits generate the highest emfs since it is here the 
tether velocity and magnetic field lines are perpendicular. 
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Other critical issues that affect the performance of an electrodynamic tether 
include plasma contactor cloud instabilities which impede the current closure process, as 
well as the characterization of the magnetosphere current closure path and its potential 
losses. The tether's susceptibility to micrometeoroid/debris damage and the effects of 
long-term insulator exposure also play a role. And finally, long tether dynamics and 
associated librations will continually affect the performance of the tether.(Bainum et al., 
1986) 
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III. THE SPACE ENVIRONMENT 



A. EARTH’S GEOMAGNETIC FIELD 

The magnetic field of the Earth is, to a first approximation, that of a sphere 
uniformly magnetized in the direction of a centered dipole axis. Both the magnitude and 
the orientation of this dipole undergo slow time variations, with a time scale of the order 
of months to years, which give rise to gradual changes of the geomagnetic field 
components. The "best fit" for this equivalent dipole model, for the smallest cumulative 
error when compared to the actual field of the Earth, is obtained by an offset of 436 km 
from the center of the Earth, displaced towards the Pacific ocean. This displaced dipole 
model, shown in Figure 8, is known as the Eccentric Dipole Model, and describes the 
measured geomagnetic field to an accuracy of two to three percent. The geomagnetic 
poles of the dipole, the location where the axis of the fictitious dipole magnet intersect the 
surface of the Earth, are located about 800 miles from the geographic poles, producing a 
tilt angle with respect to the Earth's rotation axis of 1 1.3 degrees. The north geomagnetic 
pole is located near Thule, Greenland at 78. 3N 291.0E. The geomagnetic south pole is 
located near Volstok Station, Antarctica, at 78. 3S 1 1 IE. One may accordingly define a 
geomagnetic equator as the great circle 90 degrees away from either pole in geomagnetic 
latitude. (Heinz and Olsen, 1993) 

The relationships between the geomagnetic latitude and longitude (k, A), and the 
corresponding geographic coordinates (0 0 , 4> 0 ) are given by the equations, 
sin X = cos78. 3° cos0 o cos(4) 0 - 29 1 0 ) + sin 78. 3° sin 0 O 

_ sin 78. 3° cos0 o cos(4> 0 -291°)- cos78. 3° sin 0 Q 

cos?i 
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A static dipole field can be defined by three components. With the radial, 
tangential and azimuthal components of the dipole field the magnitude and direction of the 
geomagnetic field vector can be obtained in spherical coordinates. The radial component 
of the field is given by, 

2 B 0 cos0 

B _ — 

(r/R B f 

The field's tangential component is given by, 

_ B 0 sin© 

° = (r/R c f 

The component of the field in the azimuthal direction is, 

fi =0 

Figure 9 depicts the coordinate system for the magnetic dipole field, where r is the 
distance from the dipole to the fieldpoint, R E is the radius of the Earth, 9 is the co- 
latitude of the fieldpoint, 4> is the azimuthal angle,- and B 0 is the surface field at the 



geomagnetic equator. (Tascione, 1988) 

An equation for an individual field line may be obtained by setting, 

dr _ B r 

rd& B e 

which has the solution, 

r = R 0 sin 2 0 

where R^ is the radial distance of the equatorial crossing of the magnetic field line. This 
equation, its geometry depicted in Figure 10, may be rewritten in terms of geomagnetic 
latitude, the compliment of co-latitude, 

r=R 0 cos 2 X 

The inclination of the field line, called the magnetic dip angle, is related to magnetic 
latitude by the equation, 

/ = arctan(2 tan X) 
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Finally, the magnitude of the field at any given point is obtained from solving the 
equation, 

B = — - (3sin : X + 1 y^ 1 

(r/R E y 

B. IONOSPHERIC PLASMA 

More than 99 percent of matter in the universe exists in the plasma state, a mixture 
of ionized and electrically neutral gases. A gas can be both ionized and electrically neutral 
at the same time, provided there are as many free electrons in the gas as there are net 
positive charges on the positive gas ions. Thus, plasma consists of a homogeneous 
mixture of electrons and positive ions, surging and swirling in electromagnetic bondage. 
This plasma becomes important in the Earth’s environment in the ionosphere, which 
extends upward from about fifty kilometers above the ground to about 1000 kilometers in 
altitude, and it is sustained by the ionizing action of solar ultraviolet and x-ray radiation on 
the neutral atmospheric gas. Between these altitudes, ions and electrons are present in 
sufficient quantities to affect the propagation of radio waves at frequencies from a few 
Hertz up to several megahertz.(Heinz and Olsen, 1993) 

The variation of electron density with altitude led to the subdivision of the 
ionosphere into distinctive layers - the D, E and F regions shown in Figure 1 1 . 
Ionospheric plasma is very tenuous, and is not predominantly composed of charged 
particles. Even at the most dense layers of the ionosphere, there are rarely more than one 
million electron-ion pairs in a cubic centimeter (cc) of space. In comparison, the neutral 
gas density for the same region is typically one billion particles per cc. Nevertheless, the 
positive ions and electrons are present in sufficient numbers in the ionosphere to reflect 
radio waves back to earth, and it was through this process that ground-based radio wave 
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observations first detected the region of plasma above the Earth's electrically neutral 
atmosphere. (NAS A Report, 1993) 

On the surface of the earth and upward for approximately fifty kilometers, the 
atmosphere can be considered to be an electrical insulator, and therefore there are no 
currents flowing in that region. However, convective movements of the conducting upper 
atmosphere across the earth's magnetic field lines produce electric currents in the 
ionosphere. In addition to natural variations in this plasma environment, satellites can 
alter the space plasma as they pass through their orbits. An important scientific objective 
of the study of tethers in space, involves measuring the plasma and magnetic field 
environment of an electrodynamic tether. VHF radar observations of the space 
environment can be made in an effort to characterize the environment before and after the 
satellite passes through the radar beam. These observations can be used to identify 
interactions between the tethered electrodynamic satellite system and the plasma. 

Some special features are likely to develop in the vicinity of a tether as it actively 
perturbs the plasma. Traveling at orbital velocities, the satellite affects the density, 
temperature, and electrical properties of the surrounding plasma. A plasma sheath 
develops around the satellite, while a wake, in some respects resembling a wake left by a 
boat in the water, trails behind. Due to long range interactions between charges, a plasma 
can sustain cooperative phenomena such as waves or oscillations down to much lower 
densities, and over far greater distances than ordinary gases. The plasma sheath around 
the satellite is predicted to be unstable and to change in size and shape with variations in 
ionospheric density, magnetic field alignment, and the voltage developed across the tether. 
These changes affect the type and amplitude of waves excited in the wake, creating 
observable turbulent disturbances in the environment around the satellite. (NASA Report, 
1993) 
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IV. THE PMG EXPERIMENT 



A. MISSION OBJECTIVES 

NASA's Plasma Motor Generator (PMG) mission was an active space experiment 
conducted in Low Earth Orbit (LEO), designed to validate theoretical predictions that a 
plasma discharge from an onboard hollow cathode assembly (HCA) would provide a low 
impedance electrical coupling between an orbiting spacecraft and the ambient ionosphere 
for bipolar charge transfer. In addition to the demonstration of plasma contactor 
performance was an investigation of electrodynamic-tether behavior functioning as either 
an orbit-boosting electrical motor, or as a generator that converts orbital energy into 
electricity. (Jost and Stanley, 1994) 

A comprehensive set of ground-based measurements, utilizing several 
globally-distributed ground stations, was defined to provide an evaluation of the 
interactions between the operational PMG system and the LEO environment. Sixteen 
fixed and transportable radar, magnetometer, and optical systems were used for remote 
measurements in support of the mission science objectives. This multi-sensor approach 
was used to study large-scale geophysical interactions between the orbiter and the ambient 
environment, monitor spacecraft/tether dynamics, and provide general orbital tracking 
records. Analysis of these remote measurements provided details on the following mission 

performance objectives (Jost and Stanley, 1994): 

• Far End Package (FEP) deployment to greater than 200 meters 

• Tether deployment dynamics and damping 

• Space-based hollow cathode plasma production operation 

• Current flow along geomagnetic field lines, and closure through the ionosphere 

• Tether deflection by IxB drag and thrust forces 
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B. PHYSICAL DESCRIPTION 



The experiment consisted of a tethered system of two identical plasma contactors 
connected via an 18-AWG 500 meter conducting tether. Experiment flight hardware 
included four major subsystems; the Far End Package (FEP), Near End Package (NEP), 
an electronics box, and a Plasma Diagnostics Package (PDP). The two-body tethered 
satellite system was stabilized in a gravity-gradient configuration, with the FEP oriented 
away from earth following its spring-ejection from the second stage of the DELTA-II. 
The FEP housed a hollow cathode-based plasma contactor inside an open metal box of 
dimensions 0.3m x 0.3m x 0.3m. The NEP, consisting of the electronics box and another 
plasma contactor, remained fixed to the second stage of the Delta II. The PDP was also 
attached to the rocket body, and consisted of two detectors, the Small Electron 
Electrostatic Analyzer (SESA) and the Ion Mass and Energy Analyzer (MESA).(Lilley et 
al., 1994) 

In addition to the vxB induced potential difference along the tether, batteries 
were used to bias the NEP and FEP with respect to each other. Varying resistive loads, in 
series with the tether, were cycled through during the experiment. The resulting current 
through the tether was the net effect of both contactor plasma clouds, one collecting and 
the other emitting electron current to the local space plasma environment. When 
electrons were collected at the FEP, the PMG system was in the Motor mode. (Figure 7) 
When electrons are collected at the NEP and by the Delta II rocket body, the system was 
in the Generator mode. (Figure 6) The limiting factor of this electrical circuit is the 
ionosphere’s ability to carry current. For this reason, plasma contactors must effectively 
spread the currents over a large enough area to reduce the current densities to the 
necessary levels. Thus, plasma contactors at each end of the electrodynamic tether 
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collection and emission by neutralizing space charge and scattering electrons across the 
geomagnetic field lines.(Jost and Stanley, 1994) 

Plasma contactors are a promising technology with applications in electrodynamic 
tether systems as well as grounding spacecraft to the space plasma environment. A hollow 
cathode based plasma contactor is a device that emits a dense, low-temperature plasma 
cloud through which ions and electrons are emitted and electrons are collected from the 
surrounding plasma. The construction of the hollow cathode assembly (HCA) used aboard 
PMG is illustrated in Figure 12. Functionally, the HCA heats a flow of xenon gas to 
relatively high temperatures (typically 1100° C) within a hollow electron-emitter 
(barium-oxide impregnated tungsten) cathode in the presence of a strong voltage gradient 
between the cathode and a corresponding anode plate. In this condition, the partially 
ionized gas establishes a highly-ionized "plasma-discharge," which allows current to flow 
freely between cathode element and an external anode, in this case the LEO 
ionosphere.(Jost and Stanley, 1994) 

In electron emission mode, the electrons from the plasma contactor carry the 
current while the ions neutralize the electron's space charge. In the dense contactor 
plasma, the electron density is approximately equal to the ion density. Higher energy 
electrons stream through the slowly expanding ion cloud while the lower energy electrons 
are trapped within the cloud by the potential distribution near the cathode. Each outgoing 
ion allows a number of electrons to be emitted, at roughly a ratio which is, 




When the plasma contactor is collecting electrons, the contactor ions once again 

neutralize the space charge of the ambient electrons. The contactor plasma is extremely 

k 

turbulent due to current driven electrostatic instabilities. Incoming electrons are scattered 
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by the contactor plasma and can be collected across the geomagnetic field lines within an 
effective radius of several meters. (Lilley et al., 1994) 

C. FLIGHT PROFILE 

The experiment duration, in terms of plasma contactor operation and 
consequential active environment interaction, extended six to seven hours (approximately 
four orbits) until the NEP and FEP batteries expired. As a result of the successful 
activation and extended plasma contactor operations, the ground sites typically had 
opportunities to collect data during at least three passes. Figures 13 and 14 show the 
mission flight profile and the relative positions of the four ionospheric radars used to 
view the PMG geomagnetic flux-tubes. (Jost and Stanley, 1994) 

D. RADAR OBSERVATION SITES 

The VHF radar technique utilized for PMG has been demonstrated during several 
spaceborne particle beam and plasma source experiments including NASA CHARGE-II, 
SDIO SPEAR 1, NASA ATLAS, and the Soviet-French ARAKS series. Beam steering 
monostatic (MST) doppler radars measure the doppler shift of oblique backscattered 
echoes resulting from the ubiquitous small-scale turbulence in the atmosphere. These 
phase coherent radars allow measurement of the amplitude and doppler velocity in the 
direction of probing of radio waves that are scattered back to the receiving antennas. 

Explaining MST radar observations requires two basic echoing mechanisms, 
namely scattering and reflection. For monostatic radars, backscatter and reflection arise 
from the variation of the refractive index, n, whose spatial scale along the axis of the 
radar beam is half the radar wavelength. For a near vertical beam at fifty megahertz, 
observed echoes are usually a combination of Bragg scatter (also called turbulent scatter). 
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Fresnel scatter, and Fresnel reflection. The major process causing the echoes of MST 
radars from the ionosphere is the scattering and reflection from refractive index 
irregularities due to temperature, humidity and electron density fluctuations. In addition, 
observable radar echoes from meteor trails and incoherent scatter from free electrons in 
the ionosphere should frequently occur.(Liu and Edwards, 1989) 

When a radio wave traveling at an angle to the magnetic field encounters a field- 
aligned disturbance, it is scattered such that its angle of incidence equals its angle of 
reflection. For the case when the wave is normal to the disturbance, the angle between 
the incident and reflected waves is 180 degrees and the wave backscatters. Thus a radar 
pointed such that it has a line of sight that is normal to the field lines will detect a 
maximum strength backscatter signal from any waves propagating along the field lines. 
The refractive bending of radar waves in the VHF frequencies is small enough that the ray 
paths can be approximated by straight lines.(Liu and Edwards, 1989) 

Observations of the spacecraft/environment interactions were to be obtained from 
a ground-based sensor network during the PMG mission. Scattering from refractive 
index variations caused by ionospheric plasma irregularities was expected. The 
interaction of the plasma clouds, produced by the HCA's, with the ionosphere was 
expected to cause turbulent ambient-plasma electron distributions detectable with 
ground-based VHF radar. Similarly, the refractive index of the plasma surrounding the 
operating PMG plasma source should have been overdense relative to incoming VHF 
electromagnetic waves, causing the apparent vehicle radar cross section to be evident in 
the collected radar data. Thus, two primary measurable characteristics of the PMG 
plasma clouds were expected: 

• Enhanced ionization levels and plasma turbulence in the vicinity of the vehicle 
during contactor operation. 
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• Propagation of plasma turbulence along geomagnetic flux lines over long 
distances (100's of kilometers). (Jost and Stanley, 1994) 

Four medium and large power-aperture VHF radar systems were utilized in 
support of the PMG mission, located in relative positions for favorable viewing of the 
geomagnetic flux-tubes for the first four, active-orbit ground tracks of the mission. All 
four systems operated in a pulsed, narrow-band, fixed-pointing configuration with a 
center frequency of approximately fifty megahertz. The powerful ionospheric radar at 
Jicamarca Radio Observatory in Lima, Peru (11.948 S, 76.872 W) was activated to 
observe spacecraft/environment interactions, and provide data on plasma contactor 
operation, plasma cloud size and diffusion properties, and geomagnetic field propagation 
effects. NASA/JSC deployed a transportable radar to the island of Hawaii (19.514 N, 
151.849 W) for direct observations of PMG plasma clouds, providing data 
complementary to Jicamarca. In addition, two wind-profilers, the NOAA VHF Radar 
Facility in Piura, Peru (5.167 S, 80.617 W) and the NOAA VHF Radar Facility on 
Pohnpei Island, Micronesia (6.96 N, 158.19 E) were configured for ionospheric 
measurements, supporting the Hawaii measurements, and providing a measure of 
turbulence propagated along the geomagnetic field. Finally, the USSPACECOM Radar 
Tracking network and the Kwajalein Missile Range provided precision orbital tracking of 
the tethered system during the course of the experiment. (Jost and Stanley, 1994) 

The primary goal for all the VHF radar sites was observation of 
plasma-perturbation signatures caused by PMG moving through the lower ionosphere. 
Based on the high scattering cross sections of field aligned ionospheric plasma 
turbulence, for electromagnetic radiation in the VHF bands, it was possible to observe the 
space-environmental effects of field-aligned disturbances projected great distances along 
the geomagnetic field lines. Observation of such disturbances supports the concept of 
coupled electrical current, flowing from within the conducting tether wire into the 
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rarefied ionospheric plasma. These traveling waves are theorized to stay tightly confined 
to the geomagnetic field lines until the ambient density and collision frequencies are high 
enough (at low altitudes) to allow cross-field diffusion and closure of the electrical 
circuit.(Jost and Stanley, 1994) 
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V. EXPERIMENTAL OBSERVATIONS AND ANALYSIS 



A. HAWAII PASSES 

The Transportable Radar System (TRS) deployed to Hawaii for ionospheric 
observations was a medium power-aperture coherent detection system consisting of 
several transmitter modules, two coherent multi-channel receivers, and several phased- 
array antenna systems that allowed tailoring for particular radiation pattern configurations. 
Electronic beam steering was incorporated into the design which enabled the radar to 
switch rapidly between direct illumination of the predicted PMG position in its orbit, and 
specular viewing of the projected geomagnetic lines which mapped from the orbiter. In 
this way, antenna pointing was optimized for each observation period to investigate the 
characteristics of field-aligned plasma turbulence propagating away from the active PMG 
tether-system along the geomagnetic flux tubes. 

The monostatic RX/TX phased array was capable of generating a peak power of 
160 kilowatts with a duty cycle of one percent at its operating frequency of fifty 
megahertz. An additional monopulse RX array was configured to provide interferometry 
measurements. The antenna configuration provided a gain of twenty-six decibels and was 
aimed in azimuth at the geomagnetic north pole, which for the Hawaii location was eleven 
degrees east of geographic north. It was phased in the vertical direction to produce the 
required elevation angle for the planned ionospheric observations. In its first pass over 
Hawaii, PMG was expected to pass through the radar's field-of-view at an elevation angle 
of thirty-seven degrees. Subsequent observations perpendicular to the magnetic field lines 
over Hawaii were made at an elevation angle of fifty-five degrees. (Figure 15) The radar's 
vertical and horizontal beamwidths were twenty degrees and five degrees respectfully. 
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Electromagnetic energy with a pulse width of fifty microseconds was sent out 
every interpulse period (IPP). For the first pass of PMG over Hawaii, an 1PP of five 
milliseconds was used. Following each transmission, the receivers were used to detect any 
incoming fifty megahertz radiation. Sampling every forty microseconds resulted in a range 
gate interval of 6000 meters. Due to the long slant ranges to the satellite, a coded 
waveform would be required to negate the effects of range-aliasing, while ensuring 
adequate doppler coverage, and sufficient average radiated power. However, a coded 
pulse was not used, and thus returns from two six kilometer intervals, separated by 750 
kilometers, are actually contained in each sampled gate. A total of 48 range gates were 
sampled, beginning at a slant range from the radar of 107 kilometers for PMG's first pass 
over Hawaii. 

At the Hawaii radar site, the in-phase and quadrature components of the received 
radar signal from each scattering volume were recorded continuously on analog tape. The 
vast majority of what the receivers detected was white noise background. However, if any 
turbulent propagations were present in the observed region of the ionosphere, such as 
those caused by meteor trails, then an enhancement of the incoming fifty megahertz 
radiation could be detected at a specified time delay after the initial pulse. The line of 
sight range from the radar to the returned signal could then be determined by the time 
delay between the transmitted and received signal. 

The Hawaii data underwent the first phase of processing at Systems Planning 
Corporation - Center for Space Physics in Houston, Texas. The analog tape was digitally 
sampled and analyzed, to allow plots of signal intensity as a function of time and range to 
be prepared. The power representing the signal was determined by summing the squares 
of the in-phase and quadrature components of the received signal. In order to normalize 
the received power to the system noise, the signal plus noise power was first calculated. 
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then the noise power was subtracted from it. Noise estimates were calculated by 
averaging sampled ranges that did not contain radar echoes, while signal plus noise 
estimates were calculated by averaging sampled ranges that included returns. Data 
collected during the first pass of PMG over Hawaii are shown in Figure 16. 

At the Naval Postgraduate School, additional analysis of the Hawaii data was 
performed. The backscattered signals sampled from the forty-eight range gates provided 
both intensity and doppler frequency information. A fast Fourier transform was used to 
derive this spectral information for each return echo observed in the Hawaii pass. 
Interpretation of the collected data was complicated by the background noise of meteor 
events and unstable E-layer echoes, as well as the ambiguity resulting from the range- 
aliased waveform. By comparing the spectral components of each observed return, 
however, discrimination signatures to unambiguously distinguish natural events from 
PMG-induced echoes were determined. This analysis was performed by Wayne Brewster, 
as part of his Master’s Thesis. (Brewster, 1994) The PMG satellite system was apparently 
not observed at thirty-seven degrees elevation during its first pass over Hawaii. Halfway 
through the pass, at time 14:56:00 (HHMMSS), the radar beam was switched to an 
elevation angle of fifty-five degrees to look perpendicular to the magnetic field lines. 
Approximately four minutes after the anticipated PMG meridian passage, a cluster of 
strong echoes were detected at a slant range consistent with that of the ray-path specular 
point for the geomagnetic flux-tube linked to PMG. The discrimination signatures for 
these signals were consistent with those determined to be plasma.(Brewster, 1994) 
Figures 17-19 display the spectral characteristics of these signals, as well as the time and 
slant range information about the echoes. The top panel of Figure 17 shows the intensity 
of returns in range gates 1-16. The bottom panel is a fast fourier transform versus time 
plot of the radar returns for range gate 3, corresponding to slant ranges from 869 to 875 
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kilometers. Figures 18 and 19 display similar returns for range gates 4, 5, 6, and 7, with 
each six kilometer interval beginning at a slant range of 875, 881, 887, and 893 
kilometers respectfully. A disturbance, interpreted as a traveling plasma wave can be 
identified in these figures. First seen in range gates 3 and 4, the event occurs at 
approximately time 14:58:26. Six seconds later, in range gates 5, 6, and 7, a disturbance 
is observed at time 14:58:32. 

The Kwajalein Missile Range radar was tracking the orbiting tether system as it 
passed overhead, recording the slant range, azimuth and elevation of the Delta II in its 
orbit. From these measurements an accurate orbit was derived, giving the PMG altitude 
and latitude and longitude as it orbited the earth. Projecting the orbit of PMG experiment 
forward until it reached Hawaii allowed a precise determination of the time PMG crossed 
the geomagnetic meridian to be made. Table 1 contains a portion of the orbital data for 
PMG as it made its first pass over Hawaii. Meridian crossing occurred when the azimuth 
angle reached eleven degrees, at time 14:54:19 (HHMMSS). At this time PMG was 
orbiting over the Earth at a geographic latitude of 25.8 degrees north, geographic 
longitude of 206.5 degrees east, altitude of 641.2 kilometers, and a slant range from the 
radar site of 982.5 kilometers. PMG was approaching its highest inclination on its ground 
track at this time, and was decreasing in altitude. 

While no direct scatter from the satellite system was observed at the anticipated 
meridian crossing point during this pass, through reconstruction of tracking radar data, 
the exact time PMG moved through the center of the beam and crossed the magnetic 
meridian was established. If the cluster of plasma echoes, subsequently detected with the 
radar beam perpendicular to the magnetic field, were a result of PMG’s passage through 
the field line, then these signals equate to a turbulent plasma waveform aligned with the 
magnetic field line intersected by the PMG orbiter, traveling at some propagation 
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velocity. The velocity of the plasma waveform, in the direction along the magnetic field 
line, can be calculated using the approximate flux-tube path length from PMG to the 
radar field-of view perpendicular to the field lines, and the time delay between meridian 
crossing and the detection of the plasma events. 

Using the geographic coordinates of the Hawaii radar site, its magnetic latitude 
can be calculated to be 20.5 degrees north. Likewise, using the geographic coordinates of 
the nadir point associated with PMG’s meridian crossing, a magnetic latitude of 26.3 
degrees north is found. The magnetic dip angle of the field lines at these locations can 
then be calculated using these magnetic latitudes. Above Hawaii, the field lines are 
inclined at an angle of approximately thirty-seven degrees, while at the point of meridian 
crossing this inclination increases to approximately forty-five degrees. Thus, at the point 
in the ionosphere where the Hawaii radar was pointing perpendicular to the field line, the 
value of the inclination is between these two angles. Since the inclination of the field line 
is relatively constant over this region, a straight line approximation will be used for the 
flux-tube distance between PMG meridian crossing and the field-of-view perpendicular to 
the field line. 

Forming a right triangle from the pass geometry, shown in Figure 20, the 
approximate propagation distance of the plasma waveform can be calculated by 
multiplying the slant range to PMG meridian crossing by the sine of the angle between 
the two radar beam positions. The result is a propagation distance of 303.6 kilometers. 
The time delay between meridian crossing and detection of the first plasma echo was four 
minutes and seven seconds. Using this time and distance, a propagation speed of 1 .20 
kilometers per second can be calculated. 

The nominal sources of error in this calculation include the straight line 
approximation for the geometry of the field line, six second uncertainty in time of 
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observation of the disturbance, as well as single source projection errors for the orbital 
position of the Delta II rocket based on radar tracking by Kwajalein. 

B. JICAMARCA PASSES 

In the foothills of the Andes mountains, the Radio Observatorio de Jicamarca was 
also employed to measure the intensities of field-aligned propagations. Since Jicamarca 
is located approximately one degree north of the geomagnetic equator, the local magnetic 
field lines are near horizontal, presenting an excellent opportunity to observe plasma 
propagations associated with PMG charge neutralization. With a peak power of two 
megawatts and a duty cycle of 0.6 percent, this fifty megahertz observatory is the world’s 
largest power-aperture radar. The radar’s antenna has a two-way beamwidth of one 
degree, and provides forty-five decibels of gain. A waveform with a pulse width of 6.67 
microseconds and a pulse period of 1100 microseconds was used for the ionospheric 
observations. Sampling of return echoes was divided into 250 range gates, each with an 
interval of 500 meters, providing a total window of 125 kilometers for data collection. 

Throughout the data acquisition, the Jicamarca radar was pointed perpendicular to 
the geomagnetic field lines in the lower ionosphere. The real and quadrature outputs of 
the receivers were digitized and recorded to 8-mm recording tape at the radar site. While 
this data was not available to the researchers at the Naval Postgraduate School for further 
processing, the times and ranges of initial A-scope detections of propagating plasma 
turbulence, similar to those observed in the Hawaii, were recorded and presented in the 
PMG Mission Report. (Jost and Stanley, 1994) Discrete radar detections were obtained in 
passes 2, 3, and 4, where strong echo returns from flux-tube heights consistent with PMG 
were observed. 
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Pass 2 



Pass 2 over Jicamarca was a descending pass north of the observatory (i.e., 
the satellite was moving from northwest to southeast in its groundtrack), with PMG 
losing altitude in its orbit. The tethered satellite system crossed the Jicamarca magnetic 
meridian when its reached an azimuth angle of six degrees east of geographic north. 
From the orbital tracking data in Table 2, this occurred at time 16:57:03 (HHMMSS), 
when PMG was orbiting over the Earth at a geographic latitude of 7.2 degrees south, 
geographic longitude of 283.6 degrees east, altitude of 193.2 kilometers, and a slant range 
from the radar site of 566.7 kilometers. The radar event associated with plasma 
turbulence projected along the field line that coupled the PMG system and the radar’s 
field-of-view occurred at a time of approximately 17:03:30 (plus or minus thirty seconds), 
at a slant range of 265.0 kilometers. The geometry associated with these observations is 
shown is Figure 21. 

A slightly different approach in calculating the distance traveled by the 
propagating plasma wave is used for the Jicamarca passes. Due to the dipole-behavior of 
the flux-tube geometry, the height above the radar of the field-aligned current flow is 
dependent on the latitude displacement of PMG from Jicamarca as the tethered satellite 
system crosses the magnetic meridian. Figure 22 graphs this flux-tube height dependence 
for a 300 kilometer typical-altitude reference. The height of the field line centered in the 
radar’s beam is found by adding the flux-tube height variation from Figure 22 to the 
altitude of PMG at meridian crossing. 

PMG crossed Jicamarca’s magnetic meridian during pass 2 at a time of 16:57:03. 
From the geographic position of the nadir point, the difference in latitude between 
meridian crossing and radar location was 4.7 degrees. Figure 22 indicates a flux tube 
height variation of 365 kilometers using this value. Subtracting out the 300 kilometer 
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reference height results in the magnetic field line above Jicamarca being sixty-five 
kilometers higher in altitude than at the point of PMG meridian crossing. The 
geomagnetic position of the radar site, located one degree north of the magnetic equator, 
allows the dip angle of the field lines above Jicamarca to be calculated at two degrees. 
This angle is also equal to the angle off vertical that the radar beam is pointing for 
perpendicular observation of the field lines overhead . 

From Figure 21, the arc length S can be calculated using the latitude difference 
between Jicamarca and the point of meridian crossing, and the radius from the center of 

the earth to the altitude at which PMG crossed the field line, 

S = (4.7")(— )(6378.1 km + 193.2 km)= 541.5 km 
180 

In order to find the propagation distance D, we must first subtract that part of the arc 
length calculated above that is swept out due to the inclination angle I of the radar beam, 

using a radius equal to the height at which PMG crossed the field line, 

(2")(— )(193.2 km) = 6.7 km 
180 

The difference between this two calculations is 534.7 kilometers, which approximates 
the distance the propagating plasma wave traveled in a horizontal direction along the 
magnetic field line. By forming a right triangle with this distance as one side, and the 
sixty-five kilometers traveled in the vertical direction as the other, the value of the 
hypotenuse describes the approximate total distance traveled by the wave, calculated to 
be 538.7 kilometers. 

The first radar return signals were detected at approximately 17:03:30, six minutes 
and twenty-seven seconds following PMG meridian crossing. Using this time span and 
the distance calculated above, an average propagation speed of 1.39 kilometers per 
second can be found. The thirty second uncertainty in the reported time of observance 
gives a range from 1.29 to 1.51 kilometers per second. 
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2 . 



Pass 3 



Pass 3 over Jicamarca was also a descending pass, however, this time the 
satellite system crossed the magnetic meridian south of the observatory, and was gaining 
altitude in its orbit. Meridian crossing occurred 180 degrees from the Pass 2 location, at 
an azimuth angle of 186 degrees east of geographic north. From Table 3, this occurred at 
time 18:38:29 (HHMMSS), when PMG was orbiting over the Earth at a geographic 
latitude of 17.3 degrees south, geographic longitude of 282.5 degrees east, altitude of 
215.0 kilometers, and a slant range from the radar site of 644.1 kilometers. The radar 
event associated with this pass occurred at a time of approximately 18:48:30 (plus or 
minus thirty seconds), at a slant range of 240.0 kilometers. Figure 23 illustrates the 
geometry associated with these observations. 

For Pass 3 calculations, the same method utilized in Pass 2 was used. Observe 
that PMG crossed the magnetic meridian at time 18:38:29. This time, however, PMG 
passed south of Jicamarca, and from the geographical position of the nadir point, the 
difference in latitude from the radar site was -5.4 degrees. From Figure 22, this angular 
difference equates to a flux-tube height variation of twenty-seven kilometers once the 300 
kilometer reference altitude is subtracted. Thus the field line rises twenty-seven 
kilometers from the point at which PMG crosses it to the point located at the center of the 
radar beam. 

With this pass geometry, shown in Figure 23, the desired propagation distance in 
the horizontal direction is the combination of the arc length S, and the arc length swept 

out by the inclination angle of the radar beam from local vertical, 

(5.4")(— )(6378.1 km + 215.0 km) + (2")(— )(215.0 km) = 625.4 km 
180 180 
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Once again forming the right triangle with sides of 625.4 kilometers and twenty-seven 
kilometers, the hypotenuse describes the approximate total distance of wave propagation. 
This value calculates to be 626.0 kilometers. 

The plasma echoes were reported at time 1 8:48:30 during this pass, which was ten 
minutes and one seconds following PMG meridian crossing. This time delay equates to 
an average propagation speed for this pass of 1.04 kilometers per second. The uncertainty 
again comes from the reported time of echo activity, resulting in a probable range of 0.99 
to 1.10 kilometers per second. 

3. Pass 4 

Like Pass 3, Pass 4 was also a descending pass south of the observatory, 
with PMG gaining altitude in its orbit. From Table 4, PMG reached an azimuth angle of 
186 degrees at time 20:19:1 1 (HHMMSS), marking its magnetic meridian crossing. At 
this time PMG was orbiting over the Earth at a geographic latitude of 23.7 degrees south, 
geographic longitude of 281 .9 degrees east, altitude of 282.8 kilometers, and a slant range 
from the radar site of 1365.1 kilometers. As with the previous passes, strong radar 
returns were observed, occurring at a time of approximately 20:42:30 (plus or minus 
thirty seconds), at a slant range of 445.0 kilometers. Figure 23 again illustrates the 
geometry associated with these observations. 

Pass 4 calculations are very similar to Pass 3, however, PMG crossed the meridian 
much further south in this case. At time 20:19:1 1, orbital tracking data shows that PMG 
crossed Jicamarca’s magnetic meridian. The geographic latitude associated with the 
nadir point of this crossing differed from the radar position by -1 1.8 degrees. Entering 
Figure 22 with this value results in a flux-tube height variation of approximately 181 
kilometers, once the 300 kilometer reference height is subtracted. This is the vertical 
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distance the plasma waveform transited from the meridian crossing point before passing 
through the radar beam. 

The horizontal propagation distance is calculated as for Pass 3, the combination of 
two arc lengths, 

(1 1.8" )(— >(6378.1 km + 282.8 km) + (2")(— )(282.8 km) = 1378.2 km 
180 180 

Thus, the right triangle for this pass geometry has sides of 1378.2 kilometers and 181 
kilometers, resulting in a total propagation distance of approximately 1390.0 kilometers. 

The radar event associated with this pass occurred at time 20:42:30, twenty-three 
minutes and nineteen seconds following meridian crossing. This equates to an average 
propagation speed of 0.99 kilometers per second. The thirty second uncertainty in 
detection of the event results in a probable range of 0.97 to 1.02 kilometers per second. 
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VI. SUMMARY AND DISCUSSION 



The two-body tethered satellite system that comprised NASA’s PMG experiment 
was stabilized in a gravity-gradient configuration, with hollow cathode plasma sources at 
either end of a 500 meter conducting wire. These hollow cathode plasma contactors 
coupled electric currents from the wire directly into the ambient ionospheric plasma, 
providing an excellent opportunity to observe the ionospheric modification effects, 
turbulent signatures, and spacecraft-environment interactions of an orbiting plasma 
generator. 

At any given time, the hollow cathodes at either end of the tethered satellite 
system apply potential pulses of opposite polarization to the magnetic flux tubes they 
intercept. The perturbation that occurs as a result of this charge transfer generates plasma 
waves which travel away from the PMG system, transporting the space charge set up 
between the tether end flux tubes. These turbulent plasma waveforms, triggered by the 
conducting tether and plasma contactors traversing the earth’s geomagnetic field lines in 
the ionosphere, appear to travel tightly confined along the field lines over very large 
distances. Due to the excellent scattering properties of ionospheric plasma irregularities, 
these traveling waveforms are detectable from ground-based radar. 

Through range probing of the geomagnetic flux-tubes above several VHF radar 
sites, field-aligned disturbances of a propagating plasma structure associated with PMG 
was observed on several of the spacecraft’s passes. Strong, persistent echo returns from 
the predicted flux-tube heights were observed, which mapped back to PMG at the time of 
meridian crossing with an average velocity corresponding to average ion velocities in the 
lower ionosphere. The following table summarizes the results obtained during the PMG 
experiment, from observations made in Hawaii and Jicamarca, 
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Site 


Pass 


Meridian 

Crossing 

(HHMMSS) 


Radar 

Event 

(HHMMSS) 


Time 

Delay 

(SEC) 


Propagation 

Distance 

(KM) 


Effective 

Speed 

(KM/SEC) 


HI 


1 


14:54:19 


14:58:29 


247 


303.6 


1.20 


JIC 


2 


16:57:03 


17:03:30 


387 


538.7 


1.39 


JIC 


3 


18:38:29 


18:48:30 


601 


626.0 


1.04 


JIC 


4 


20:19:11 


20:42:30 


1399 


1390.0 


0.99 



The observation of two distinct signals during the Hawaii pass remains an unsolved 
problem. Perhaps a small component of the wave's velocity in the radial direction causes 
the echoes to be observed at the two specified times. Notice also that some diffusion or 
eddying motion occurs over the large distances of observation involved in this 
experiment, seen by the decreasing waveform propagation speeds in the successive 
Jicamarca passes as the propagation distances increased. 

In comparison to the average propagation velocities shown in the table above, the 
thermal velocity of the ionosphere electrons is given by the equation. 




where k is Boltzmann’s constant (1.38x 10' 23 J /" K ), T e is the average temperature of the 
electrons (use 1000°K), and m e is the electron mass (9.1x10 }] kg). Calculating this 
value results in an electron thermal velocity of approximately 123 kilometers per second - 
much greater than the observed velocities above. Assuming the thermal temperature of 
the ions is the same as the electrons, and substituting the oxygen ion atomic weight 
(15.999) into the above equation, the ion thermal velocity can be calculated to be 
approximately 0.72 kilometers per second. This value is of the same order as the 
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propagation velocities observed in the table above. Since the hollow cathode plasma 
sources were emitting xenon gas, this calculation is repeated using xenon ions (atomic 
weight of 131.30). The resulting xenon ion thermal velocity is 0.25 kilometers per 
second, which is close to an order of magnitude less than the values in the table. 

Plasma waves propagating away from an electrodynamic tether system in low- 
earth orbit have been postulated to be Alfven and fast-mode waves. Alfven wave 
dispersion can be thought of as representing the wave that propagates when a field line, 
loaded with plasma, is plucked in the transverse direction. Such waves travel tightly 
confined to the Earth’s magnetic field lines, and possible reflection of these waves may 
occur at density gradients like the E-layer of the Earth’s ionosphere. (Dobrowolny, 1993) 
The Alfven velocity is given by the equation, 

V A = -‘—y 

where fj. () is the permeability constant (4nx 10“ 7 H / m), p m = is the ion mass density, 
and the measured value of the Earth’s magnetic field at the Earth’s surface at the equator 
is *.=3.13 x 10 5 7\ Using a value for the ionospheric number density of n n = \0" m 3 

and assuming oxygen ions, the resulting Alfven velocity is 540 kilometers per second. 
Repeating the calculation for xenon ions results in an Alfven velocity of 189 kilometers 
per second. Neither of these results are of the same order as the observed propagation 
velocities above Hawaii and Jicamarca. 

Following the results of Urrutia et al., 1994, it may be that these plasma 
disturbances propagate in a whistler mode. These results differ fundamentally from the 
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traditional current model: instead of Alfven wings and a phantom current loop, Urrutia 
predicts whistler wings and a short diffuse current loop, due to cross-field shunting 
electron Hall currents. At each position in their orbit, the plasma contactors on the 
electrodynamic tether excite a low-frequency whistler wave packet along the magnetic 
field line. The superposition of these wave packets forms a coherent, wing-like current 
structure propagating away from the tether system, which travel along the field lines at 
whistler speeds and disperse by inducing secondary plasma currents. The results of this 
study suggest that an electrodynamic tether cannot generate a long, filamentary current 
loop. Instead, the current loop is shunted by cross field currents associated with the 
continuous shedding of whistler waves, well before encountering the boundary of the 
plasma. Such a whistler wing generated by a long space tether could possibly produce a 
whistler tone detectable at the ground when the tether passes by. (Urrutia et al., 1994) 

In order to find the phase velocity of a whistler wave traveling though the 
ionosphere, it is necessary to first solve for the refractive index of the ambient plasma.. 
This refractive index, n, can be found from the following equation, 



where o) is the frequency of the wave, the plasma frequency, (o p , is given by the equation. 



and the electron gyrofrequency, co v , is given by the equation, 
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For these equations, n 0 is the ionosphere number density, e is the electron charge 
(1.6x 10" I9 C), £„ is the permittivity constant (8.85x 10' l2 F/m), m e is the electron mass, 
q is the charge on the particle - in this case an electron, and B is the strength of the 
magnetic field (use B 0 ). Calculating the plasma frequency and the electron 
gyrofrequency using these numbers results in values of 2.84 megahertz and 0.88 
megahertz, respectfully. For the special case of whistler wings, where co«co c , the 
square of the refractive index of the plasma is inversely proportional to the frequency of 
the wave, 



n 2 = 1+- 



C o„ J__ (0^ 

ft), ft) ft), ft) 



Finally, the phase velocity of the traveling wave can be calculated using the relationship. 



V=“ 

n 

where c is the speed of light in a vacuum (3.0x 10 8 m/.y). Considering possible wave 
frequencies from 0-1000 Hertz, the corresponding phase velocity was calculated and 
plotted in Figure 24. For frequencies in the ELF/VLF range, the phase velocity of 
whistler wings appear to be in the 100-1000 kilometer per second range, again far greater 
than the propagation velocities in the table above. 

In conjunction with the PMG radar measurements, the University of Genoa, Italy 
provided a highly sensitive ground-fixed magnetometer for observations at the Earth’s 



45 



surface of spontaneous and man-made emissions, propagating from the orbiting tether 
through the Earth-ionosphere cavity. This super-conducting quantum interference device 
(SQUID) magnetometer was used to investigate perturbations in the ambient natural 
magnetic fluctuations during PMG’s pass overhead. Although final analysis has not yet 
been completed, preliminary results so far do not exclude the presence of possible 
satellite perturbations, such as whistler wings, affecting the ELF background noise. 
(Minna et al., 1994) 

The table below summarizes the various propagational wave calculations made in 
the preceding paragraphs. Only the oxygen ion thermal velocity appears to be of the same 



electron thermal velocity 


123 km/sec 


0+ ion thermal velocity 


0.72 km/sec 


Xe+ ion thermal velocity 


0.25 km/sec 


0+ Alfven velocity 


540 km/sec 


Xe+ Alfven velocity 


189 km/sec 


ELF/VLF whistler wave velocity 


1 00- 1 000 km/sec 



order of magnitude as the proposed propagation velocity of the plasma disturbances 
observed at Hawaii and Jicamarca. This suggests that the traveling disturbances 
associated with PMG’s meridian crossing are ion sound waves, composed of the ambient 
plasma of the ionosphere. Alternatively, these signals might be thought of as ambipolar 
diffusion of plasma from the cathodes along the field line at approximately the ion sound 
speed. (Chen, 1984) 
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VII. CONCLUSIONS 



Multiple observations of the plasma-perturbation signatures associated with PMG 
meridian crossing supports the concept of disturbances, caused by passage of an 
electrodynamic tether through the ionospheric plasma, propagating along the magnetic 
field lines. In each instance, the observed plasma disturbance lasted several seconds and 
had a propagation velocity of the order of 1000 meters per second. The large distances of 
propagation suggest that the signals would continue their path along the field lines, until 
reaching the E layer of the ionosphere. If these signals represent field-aligned currents, 
this would allow for current loop closure. 

The true dynamics of the phantom current loop still remain a mystery, however, 
further processing of all data collected for each pass during this experiment will provide 
more detail regarding signal strength, spectral behavior, spatial diffusion, and echo 
duration. This detail should be more interpretable in terms of field-aligned current flow, 
current closure phenomenology, and general electrified-spacecraft interactions with the 
ionospheric environment. The complex sheath structure around a highly charged body 
moving in a magnetoplasma and the processes occurring there, such as wave generation 
and plasma expansion phenomena, are still largely unexplored. A large wake is probably 
characterizing the satellite environment, and wave and particle phenomena associated 
with this wake, and with the expansion of the outside plasma into this evacuated cavity 
needs to be investigated. 
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The concepts for use of tethers in space have been with us for some time. 
However, the continuing new developments concerning electrodynamic tethers now 
provide a revolutionary capability for generating propulsion through momentum exchange, 
and in generating electric power through interactions with the space environment. 
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Figure 1: Early Gravity-Gradient Concepts 



49 



/ 




Local Vertical 




Figure 2: Forces Acting on Tethered Satellites 
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Figure 3: Restoring Forces Acting on Tethered Satellites 
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Figure 4: Current Paths for Electrodynamic Tethers 
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Figure 5: The Phantom Current Loop 
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Figure 6: The Generator Principle 
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Figure 7: The Motor Principle 
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Figure 8: The Eccentric-Dipole Model of the Earth’s Magnetic Field 
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Figure 9: Geometry for a Static Dipole Field 
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Figure 10: Geomagnetic Latitude and Equitorial Crossing Distance of a Field Line 
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Figure 11: Ionospheric Electron Density Layers 
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Figure 12: PMG Hollow Cathode Assembly Schematic 
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Figure 13: PMG Mission Flight Profile 
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Figure 14: PMG Ground Track and Observing Radar Sites 
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Figure 15: Hawaii Radar Beam Geometry for Pass 1 
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Figure 16: PMG Hawaii Pass 1 Echo Activity (Jost and Stanley, 1994) 
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PMG Pass 1: Time 14 58 





50 data samples, zero padded for 256 element fft 



200 

150 

100 

50 

o 



Figure 17: Hawaii Pass 1 Plasma Returns - Range Gate 3 
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Figure 18: Hawaii Pass 1 Plasma Returns - Range Gates 4 and 5 
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Figure 19: Hawaii Pass 1 Plasma Returns - Range Gates 6 and 7 
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Figure 20: Geometry for Hawaii Pass 1 Calculations 
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Figure 21: Geometry for Jicamarca Pass 2 Calculations 
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FLUX-TUBE MAPPING OVER JICAMARCA 
AT GEOMAGNETIC MERIDIAN 




70 



Latitude Difference (degrees) 



Propagation Distance 




Arc Length 



PMG 



Field 

Line 




Figure 23: Geometry for Jicamarca Pass 3 and 4 Calculations 
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Figure 24: Formation of Alfven Wings 
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Figure 25: Formation of Whistler Wings 
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Figure 26: Whistler Mode Phase Velocity/Frequency Relationship 
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TABLE 1: PMG MERIDIAN CROSSING FOR HAWAII PASS 1 
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93177201916.920 -23.7930 282.2805 284.408 5.9534 184.0800 1369.767 

93177201917.940 -23.8054 282.3504 284.670 5.9524 183.7699 1370.679 

93177201918.960 -23.8178 282.4203 284.932 5.9508 183.4601 1371.628 

93177201919.980 -23.8302 282.4902 285.194 5.9488 183.1509 1372.614 
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